Systematic Analysis of Majorization in Quantum Algorithms 
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Motivated by the need to uncover some underlying mathematical structure of optimal quantum 
computation, we carry out a systematic analysis of a wide variety of quantum algorithms from the 
majorization theory point of view. We conclude that step-by-step majorization is found in the known 
instances of fast and efficient algorithms, namely in the quantum Fourier transform, in Grover's 
algorithm, in the hidden affine function problem, in searching by quantum adiabatic evolution and 
in deterministic quantum walks in continuous time solving a classically hard problem. On the 
other hand, the optimal quantum algorithm for parity determination, which does not provide any 
computational speed-up, does not show step-by-step majorization. Lack of both speed-up and step- 
by-step majorization is also a feature of the adiabatic quantum algorithm solving the 2-SAT "ring 
of agrees" problem. Furthermore, the quantum algorithm for the hidden affine function problem 
does not make use of any entanglement while it does obey majorization. All the above results give 
support to a step-by-step Majorization Principle necessary for optimal quantum computation. 

PACS numbers: 03.67.-a, 03.67.Lx 



I. INTRODUCTION 

One of the main open problems in quantum computa- 
tion theory is finding some mathematical structure un- 
derlying optimal quantum algorithms. There is a rather 
short list of ideas behind the design of fast algorithms. 
Grover's quantum searching algorithm pj exploits calls 
to an oracle by enhancing a particular state via a rota- 
tion in the relevant Hilbert space associated to the prob- 
lem. Shor's quantum factoring algorithm Q exploits the 
periodicity of an initial quantum state using only a min- 
imum of Hadamard gates at the core of the quantum 
Fourier transform. Based on more general quantum me- 
chanical principles, the idea of using adiabatic evolution 
to carry quantum computation Q has proven suitable 
for performing Grover's algorithm and has been numer- 
ically studied as a candidate for attacking NP-complete 
problems. More recently, deterministic random walks in 
continuous time have been proven to solve a classically 
hard problem in polynomial time Q. Many other quan- 
tum algorithms can be mapped to the above families and, 
therefore, bring no further insight. 

Some attempts to uncover an underlying principle, 
common to all known optimal algorithms, have already 
been explored though not definite and satisfactory an- 
swer has been found. One relevant instance is the role of 
entanglement in quantum algorithms 0, El IS HI El • 
Although entanglement is a natural resource to be ex- 
ploited in quantum algorithm design, there are known 
examples of fast algorithms where the quantum register 
remains in a product state all along the computation. 
Our work will concentrate on quite a different proposal. 
The basic idea is that efficiency is related to a strong 
majorization principle. We shall investigate the way the 
probability distribution associated to the computational 
basis evolves along optimal quantum computations and 
find that it obeys a very constrained behavior we shall 
analyze in detail. 



Let us recall that majorization theory arises as the nat- 
ural framework to study the measure of disorder for clas- 
sical probability distributions [HE H, III The no- 
tion of ordering emerging from majorization is far more 
severe than the one quantified by the standard Shan- 
non entropy. If one probability distribution majorizes 
another, a set of inequalities must hold to constrain the 
former probabilities with respect to the latter. These in- 
equalities lead to entropy ordering, but the converse is 
not necessarily true. In quantum mechanics, it has been 
proven that majorization is at the heart of the solution of 
a large number of quantum information problems. Ma- 
jorization plays a fundamental role in topics like ensemble 
realization, conversion of quantum states via local opera- 
tions and classical communication, and characterization 
of positive operator valued measurements [T3 | . 

In the context of quantum algorithms, a majorization 
principle has been formulated, proven for Grover's algo- 
rithm and verified for Shor's algorithm in Ref. |l6l |. Fur- 
thermore, a complete proof of majorization in quantum 
phase-estimation algorithms was presented in [T3. The 
underlying idea behind these analysis is that the prob- 
ability distribution associated to the quantum state in 
the computational basis is step-by-step majorized until 
it is maximally ordered. Then a measurement will pro- 
vide the sought solution with high probability. It has 
also been proven that the way such a detailed majoriza- 
tion emerges in both algorithmic families is intrinsically 
different [L7I ] . 

In this paper we analyze the consistency and univer- 
sality of a possible Majorization Principle. More con- 
cretely, we have studied several different quantum algo- 
rithms based on distinct properties. First, we have con- 
sidered the problem of finding a hidden affine function 
by means of calls to an oracle. This problem is relevant 
because it is a known fast quantum algorithm that uses 
no entanglement at all. Second, we have taken a non- 
efficient instance, namely the parity determination prob- 
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lcm. This is a case where the final solution must match a 
global majorization, yet it does not obey step by step ma- 
jorization neither presents any quantum speed-up. The 
third case considered here is the large class of quantum 
adiabatic evolution algorithms. Efficiency and optimal- 
ity has been proven to depend on the interpolating time 
path taken along the evolution. It is a remarkable fact 
that optimality appears when step-by-step majorization 
is present. Finally, deterministic quantum random walks 
provide exponential speed-up over classical oracle based 
random walks. Again, a manifest strong majorization 
arrow is detected. 

The conclusion of the accumulated research is that all 
the analyzed instances of quantum algorithms support 
a step-by-step Majorization Principle. That is, opti- 
mal quantum computation is systematically verified to 
correspond to a step by step detailed reordering of the 
whole probability distribution in the computational ba- 
sis. Some of the instances show the extra feature that 
the initial state can be prepared in different ways. Then, 
an initial step-by-step reverse majorization period pre- 
cedes the subsequent step-by-step majorization, closing 
an invertible cycle. The study of quantum computation 
by adiabatic evolution shows the possibility of slower al- 
gorithms that maintain majorization. This implies that 
step-by-step majorization may be a necessary but is def- 
initely not sufficient condition for efficiency 

We have structured the paper as follows: in Sec. II we 
briefly review some concepts about majorization theory 
and how it is related to quantum algorithms. We develop 
an analysis of a quantum algorithm for solving a hidden 
affine function problem in Sec. III. In Sec. IV we study 
an optimal quantum algorithm solving the parity prob- 
lem. We move to an investigation of adiabatic quantum 
computation in Sec. V, analyzing the effect of the evolu- 
tion path in adiabatic searching algorithms in Sec. V.A 
and Sec. V.B, as well as the effect of the speed in the 
time variation of the Hamiltonian in Sec. V.C. A further 
example of adiabatic evolution solving the 2-SAT "ring of 
agrees" problem is provided in Sec. V.D. In Sec. VI we 
examine a recently proposed quantum algorithm based 
on a continuous time quantum walk on a graph solving a 
classically hard problem. Finally, in Sec. VII we state a 
Majorization Principle based on all the previous results 
and collect our conclusions. 



A. Majorization theory 

Let us consider two vectors x, y £ R d such that 

^2i=i x i = 2i=i Vi ~ 1; whose components represent 
two different probability distributions. We say that dis- 
tribution y majorizes distribution x, written as x -< y, if 
and only if there exist a set of permutation matrices Pj 
and probabilities pj such that 

3 

Because the probability distribution x can be obtained 
from y by means of a probabilistic sum, the definition 
given in equation JQ) provides the intuitive notion that 
the x distribution is more disordered than y. An alterna- 
tive and usually more practical definition of majorization 
can be stated in terms of a set of inequalities to be held 
between the two distributions. Consider the components 
of the two vectors sorted in decreasing order, written as 
(zi, . . . Zd) = vr. We say that y^ majorizes x^ if and only 
if the following relations are satisfied: 

k k 

Xi < J2 Vi k = l...d. (2) 

In this paper we call probability sums similar to the ones 
appearing in the previous expression as "cumulants" . 
There is still a third way of defining majorization involv- 
ing the use of doubly stochastic matrices. A real d x d 
matrix D = (-Dy) is said to be doubly stochastic if it 
has non- negative entries, and each row and column of D 
sums to 1. We say that y majorizes x if and only if there 
is a doubly stochastic matrix D such that 

x = Dy . (3) 

The equivalence among the three given definitions can 
be proven [T^ |. Complementarity, we say that the prob- 
ability distribution x reversely majorizes distribution y 
if and only if y majorizes x. A powerful relation in- 
volving majorization and the common Shannon entropy 
H(x) = — X)f=i x i log^i °f probability distribution x is 
that if x -t, y then H(y) > H(x). 



II. MAJORIZATION THEORY AND ITS 
RELATION TO QUANTUM ALGORITHMS 

Our approach to the mathematical study of quantum 
processes will be through majorization's eye. We review 
in this section the contact point between majorization 
theory an d q uantum algorithms, as stated previously in 
and |l7j . In particular, we also present here the 
concept of "natural majorization", first stated in [l7j . 
which will eventually be used in this work. 



B. Link between majorization theory and quantum 
algorithms 

The way we relate majorization theory to quantum 
algorithms can be stated as follows: let \tp ijn ' > ) be the pure 
state representing the register in a quantum computer at 
an operating stage labeled by m = 1 . . . M, where M 
is the total number of steps in the algorithm, and let 
N be the dimension of the Hilbert space. If we denote 
as {|i)}£Li the basis in which the final measurement is 
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performed in the algorithm, we can naturally associate a 
set of sorted probabilities p il i = 1 . . . N, to this quantum 
state in the following way: decompose the register state 
in the measurement basis such that 

|^>> = f>^|<>. (4) 

i=l 

The probability distribution associated to this state is 

P W = {p! m) } pH^|a^| 2 = |(#^)| 2 (5) 

where i = 1 . . . N, which corresponds to the probabili- 
ties of all the possible outcomes if the computation were 
stopped at stage to and a measurement were performed. 
A quantum algorithm will be said to majorize this prob- 
ability distribution between steps to and m + 1 if and 

only if mini 

p M ^ p (m+l) _ ( 6 ) 

Similarly, a quantum algorithm will be said to reversely 
majorize this probability distribution between steps m 
and m + 1 if and only if 

p (™+i) ^ p M . ( 7 ) 

If the situation given in equation JBJ) is step-by-step 
verified, there is a net flow of probability directed to the 
values of highest weight, in such a way that the probabil- 
ity distribution will be steeper as time flows. In physical 
terms, this can be stated as a very particular constructive 
interference behavior, namely, a constructive interference 
that has to step-by-step satisfy the constraints given in 
equation @. The quantum algorithm builds up the so- 
lution at each time step by means of this very precise 
reordering of probability distribution. 

It is important to note that majorization is checked on 
a particular basis. Stcp-by-step majorization is a basis 
dependent concept. Nevertheless there is a preferred ba- 
sis, which is the basis defined by the physical implemen- 
tation of the quantum computer or computational basis. 
The principle we analyze is rooted in the physical pos- 
sibility to arbitrarily stop the computation at any time 
and perform a measurement. The claim pursued along 
the paper is that the probability distribution associated 
to this physically meaningful action obeys majorization. 

C. Natural majorization in quantum algorithms 

Let us now define the concept of natural majorization 
for quantum algorithms as it was originally presented in 
[l7| . Working with probability amplitudes in the basis 
{|*)}fLi? the action of a particular unitary gate at step 



to makes the amplitudes evolve to step to + 1 in the 
following way: 

N 

^ n+1) = E M m) < (8) 

where Uij arc the matrix elements in the chosen basis of 
the unitary evolution operator (namely, the propagator 
from step m to step m + 1). Inverting the evolution, we 
can write 

N 

4 m) = E^4 m+1) > ( 9 ) 

3=1 

where Cij are the matrix elements of the inverse unitary 
evolution (which is unitary as well). Taking the square 
modulus we find 

N 

| a M| 2 = |Cy | 2 | a J" +1) | 2 + interference terms . (10) 

3=1 

Should the interference terms disappear, majorization 
would be verified in a "natural" way between steps m 
and to + 1 because the initial probability distribution 
could be obtained from the final one only by the action 
of a doubly stochastic matrix with entries |CV, | 2 . This is 
the so-called "natural majorization" : majorization which 
naturally emerges from the unitary evolution due to the 
lack of interference terms when making the square mod- 
ulus of the probability amplitudes. Similarly, we can de- 
fine the concept of "natural reverse majorization" , which 
follows in a similar way: there will be "natural reverse 
majorization" between steps to and m + 1 if and only if 
there is "natural majorization" between time steps m+1 
and m. 



D. Majorization in Grover's and Shor's quantum 
algorithms 

In order to motivate the forthcoming study we briefly 
sketch here the main results found concerning the anal- 
ysis of majorization in the two well-known quantum al- 
gorithms of Grover 0] and Shor B. These results were 
already presented in Ref. and |l7j , so we address the 
reader interested in precise details to these references. 

On the one hand, Grover's quantum searching algo- 
rithm was found in [T^ | to follow a step-by-step ma- 
jorization. More concretely, each time the Grover's op- 
erator is applied the probability distribution obtained 
from the computational basis obeys the set of constraints 
given in equation J2J until the searched state is found. 
Furthermore, because of the possibility of understand- 
ing Grover's quantum evolution as a rotation in a two- 
dimensional Hilbert space (see for example [2^') it is 
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seen that the quantum algorithm follows a stcp-by-step 
reverse majorization when evolving far away from the 
marked state, until the initial superposition of all possi- 
ble computational states is obtained again. The quantum 
algorithm behaves such that majorization is present when 
approaching to the solution, while reverse majorization 
appears when escaping from it. A cycle of majorization 
and reverse majorization emerges in the process as we 
consider long enough time evolutions, due to the rota- 
tional nature of Grover's operator. 

On the other hand, Shor's quantum algorithm was 
analyzed inside of the broad family of quantum phase- 
estimation algorithms. In 16] it was observed that a 
step-by-step majorization seemed to appear under the ac- 
tion of the last quantum Fourier transform when consid- 
ered in the usual Coppersmith decomposition j2j|. One 
step further was taken in |l7j . were the complete math- 
ematical proof of this property was provided. The result 
relies on the fact that those quantum states that can 
be mixed by a Hadamard operator coming from the de- 
composition of the quantum Fourier transform only differ 
by a phase all along the computation. Such a property 
entails as well the appearance of natural majorization, 
in the way presented above. Natural majorization was 
proven relevant for the case of Shor's quantum Fourier 
transform. This particular algorithm manages stcp-by- 
step majorization in a most efficient way. No interference 
terms spoil the majorization introduced by the natural 
diagonal terms in the unitary evolution. It is still unclear 
the role that natural majorization plays in distinguishing 
different level of efficiency in quantum algorithms. 

These two results suggest a possible relation between 
majorization and quantum algorithms. This is the point 
to be exploited in detail with further examples in the 
next sections. 



III. ANALYSIS OF A QUANTUM ALGORITHM 
FOR SOLVING A HIDDEN AFFINE FUNCTION 
PROBLEM 

The problem of finding hidden affinc functions was ini- 
tially proposed by Bernstein and Vazirani [Tsj as a gen- 
eralization of Deutsch's problem ^||- Further studies 
have investigated into this class of problems, providing 
a range of fast quantum algorithms for solving different 
generalizations |2(J, |21j . The case we present in this work 
is one of the multiple variations stated in Ref. |2l|, but 
the main results we obtain can be verified as well for the 
whole set of quantum algorithms solving similar prob- 
lems. 



x, to, 6 €E Zn, find to. 

The classical analysis reveals that no information 
about to can be obtained with only one evaluation of the 
function /. Conversely, given the unitary operator Uf 
acting in a reversible way in the Hilbert space H^v <£> Hjv 
such that 



U f \x)\y) = \x)\y + f{x)) , 



(11) 



(where the sum is to be interpreted as modulus N) , there 
is a quantum algorithm solving this problem with only 
one query to Uf. 



B. Quantum algorithm 

Let us take N = 2™, being n the number of qubits. 
The quantum algorithm optimally solving the problem 
previously presented reads as follows: 

• Prepare two registers of n qubits in the state 
|0...0}|V>i) € Hat <g> Hjv, where |^i) = 
QF 1 T(A^)- 1 |1), and QFT(N)~ 1 denotes the inverse 
quantum Fourier transform in a Hilbert space of di- 
mension N. 

• Apply QFT(N) over the first register. 

• Apply Uf over the whole quantum state. 

• Apply QFT(N)~ 1 over the first register. 

• Measure the first register and output the measured 
value. 

The different steps concerning this process are summa- 
rized in Fig. ^ 
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FIG. 1: Quantum circuit solving the hidden affine function 
problem. Each quantum wire is assumed to be composed of 
n qubits. The arrow at the end indicates a measurement. 



A. Setting of the problem 

Let us state the following problem (see pljl: 

Given an integer N function f : x — > mx + b, where 



C. Analysis of the quantum algorithm 

We now show how the proposed quantum algorithm 
leads to the solution of the problem. Our analysis raises 
two observations concerning the way both entanglement 
and majorization behave in the computational process. 
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In the first step of the algorithm, the quantum state 
is separable, noting that the quantum Fourier transform 
(and its inverse) applied on a well defined state in the 
computational basis leads to a perfectly separable state 
(see, for example, |2(j) Actually, this separability holds 
also step-by-step when a decomposition for the quan- 
tum Fourier transform is considered, such as the Copper- 
smith's decomposition [23j. That is. the quantum state 
1 ... 0) is unentangled. 

The second step of the algorithm corresponds to a 
quantum Fourier transform in the first register. This 
action leads to a step-by-step reverse majorization of the 
probability distribution of the possible outcomes while 
it does not use neither create any entanglement. More- 
over, natural reverse majorization is at work due to the 
absence of interference terms. 

Next, it is easy to verify that the quantum state 

1 £z? 

l^> = 7^E e ~ 2 ^'> ( 12 ) 

V 3=0 

is an cigenstate of the operation \y) — > \y + f{x)) with 

eigenvalue e 2 ™ - " - . After the third step, the quantum 
state reads 

cc — \ — / 

(13) 

The probability distribution of possible outcomes has not 
been modified, thus not affecting majorization. Further- 
more, the pure quantum state of the first register can be 
written as QFT(N)\m) (up to a phase factor), so this 
step has not created any entanglement among the qubits 
of the system either. 

In the fourth step of the algorithm, the action of the 
operator QFT(N)~ 1 over the first register leads to the 
state e 27 ™"« |to)|i/;i). A subsequent measurement in the 
computational basis over the first register provides the 
desired solution. Recalling the results found in 0, we 
see that the inverse quantum Fourier transform natu- 
rally majorizes step-by-step the probability distribution 
attached to the different outputs. On the other hand, the 
separability of the quantum state still holds step-by-step. 
Note that step-by-step majorization is in fact dependent 
on the specific implementation of the quantum Fourier 
transform operation, but nevertheless it holds true for 
other decompositions of the operator appart from the 
usual Coppersmith's one (as we already stated in [l7j K 

It is clear that the quantum algorithm is faster than 
any of its possible classical counterparts, as it only needs 
of a single query to the unitary operator Uf to get the 
solution. We can summarize our analysis of majorization 
for the present quantum algorithm as follows: 

Result 1: The fast quantum algorithm for finding a 
hidden affine function shows a majorization cycle based 
on the action of a QFT(N) and a QFT(N)- 1 . 



We understand that the algorithm is entanglcmcnt-free 
as long as we analyze it between the action of the differ- 
ent unitary gates. From a more physical point of view, 
the quantum registers may become highly entangled dur- 
ing the performance of multi-qubit gates, despite it is not 
present between two of them. Our assertion relics then 
on the study of the system at these particular steps in 
the computation, which wc think to be the most natural 
steps to consider. It follows that there can exist quantum 
computational speed-up without the use of entanglement 
(in the way made precise before). In this case, it is seen 
that no resource increases exponentially. Yet, a majoriza- 
tion cycle is present in the process, which is indeed rooted 
in the structure of both the quantum Fourier transform 
and the quantum state. 

IV. ANALYSIS OF AN OPTIMAL QUANTUM 
ALGORITHM FOR SOLVING THE PARITY 
PROBLEM 

The problem of finding the parity of a given function 
/ : x € Jjn — > { — 1, +1}, usually known as the parity 
problem, has been shown to be as hard for a quantum 
computer in the quantum oracular model of computa- 
tion as it is for a classical computer 0, HH : no quantum 
speed-up can be achieved in this case. We shall first 
present the problem and then analyze an optimal quan- 
tum algorithm proposed in [24[ from the point of view of 
majorization. 

A. Setting of the problem 

Let us state the parity problem in the following way 
(see 0): 

Given a function f : x G 1>n — * { — 1, evaluate the 

product of all the f{x) for all the possible x. 

It has been proved that a quantum computer will need 
at least N/2 queries to the quantum oracle / for sol ving 
this problem compared to the N classical queries [24| . 
Thus, a quantum computer is not faster than a classi- 
cal one (in the limit of a very large input, where N goes 
to infinity) when dealing with this situation: no better 
efficiencies can be obtained using the quantum computa- 
tional paradigm in getting the desired result. The time 
complexity of the best possible quantum algorithm will 
be O(N), without improvement with respect to the clas- 
sical time complexity, because the quantum speed-up is 
only by a factor of two. 

B. Quantum Algorithm 

Let us briefly outline the main lines of an optimal quan- 
tum algorithm solving the parity problem in exactly N/2 
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queries to the oracular function /, which was initially 
presented in [24||. We first introduce a series of defini- 
tions and notations: the function f{x) will be evaluated 
through a quantum oracle acting in the following way 



U f \x,+l) = \xJ(x)) 



U f \x,-1) = \x,-f(x)) , 



(14) 



where the second register is a qubit taking the values 
±1. Let us define also the quantum state \x,a) = 
(\x, +1) — \x, —1)), which is seen to be a proper state 
of operator Uf with eigenvalue f(x). With this defini- 
tions, the quantum algorithm reads as follows: 



• Prepare the initial quantum state \ipo) 
l 



1 ! Ex=i \ x ' a )- 



• Apply the following operations over the initial 
quantum state: 



V N/2 UfV N/2 -iUt ■ ■ ■ ViUffa) = |V/) , (15) 

where U / is defined as before, and the rest of oper- 
ators are defined as Vat/2 = 1, ^ = 7 V i ^ N/2, 
with 



N 
~2 



V\x,a) 


= \x + 


1, a) x 


V\f,a) 


= im: 


) 


V\x,a) 


= |a; + 


1, a) x 


V\N,a) 






= 1-- 
1 2 


f l,o) . 



N 
~2 



(16) 



• Measure the observable \ipo}(ipo\ over \ipf) . 

Note that the final measurement is to be made on 
a specific basis. 

C. Analysis of the quantum algorithm 

We now make an study of how this algorithm leads 
to the desired solution of the proposed problem. The 
analysis will show the way majorization behaves in this 
optimal but non-efficient quantum process. 

We focus on how the operations in the second step 
affect the quantum state leading to the solution, and thus 
affect the probability distribution of possible outcomes 
for the final measurement. If we apply Uf to the initial 
state, we get 



1 N 

UfWo) = -7= J £f( x )\ x > a ) 

v x — l 



After the application of operator Vi the quantum state 
evolves to 



N/2 

nufito) =- 7 =][]/(*-i)M 

V x — l 

1 N 

+7^ E (is) 

V iV x=N /2+l 

as can be directly checked (care must be taken with the 
possible values of x in both sums). If we now apply again 
Uf we get: 

N/2 

UfViUf\tM =-=^/(x)/(.T-l)|x,a) 

V X—l 

1 N 

7= E mf(x-l)\x,a) (19) 



x=N/2+l 

so we begin to recognize the pattern the algorithm fol- 
lows. At the end of the computation the final state is 



N/2 

=/(l).../(JV/2)5>,o> 

X — l 

N 



(17) 



+-^/(JV/2+l).../(JV) J2 l^ a )'( 2 °) 

and it is easily verified that this state is equal to \tpo} if 
the parity of the function is equal to +1, and orthogonal 
to jVo) (say \i/Jq )) if the parity is — 1. A suitable mea- 
surement of the observable \ipo)(tpo\ can then distinguish 
between the two values. 

Let us now analyze the way this algorithm behaves 
with respect to majorization. As stated in Sec. II. B, 
majorization must always be checked from the proba- 
bility distribution of obtaining the different outcomes of 
the final measurement. In other words, the probability 
distribution subject of analysis must always be the one 
obtained from the final measurement basis. Such a basis 
turns usually to be the computational one, but it is not 
necessarily so. Here, we are dealing with one of these 
exceptional cases in which the final measurement basis 
differs from the computational one. Consequently, ma- 
jorization must be studied in this unusual but natural 
basis for the quantum algorithm. 

The only two vectors we know of this basis are |Vo) 
and \tpQ ). We could extend them to a whole basis but it 
is not necessary for our purposes, as we can analyze the 
probability of being in each of these two states all along 
the computation. Should majorization be present step- 
by-step in the process, the probability of being in one of 
these two states would smoothly decrease in favor of the 
other one, which would parallely smoothly increase. In a 
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naive way, this is what a majorization arrow means: as 
the process evolves the probability of being in the right 
state becomes bigger and bigger. 

We can observe that this does not happen in the al- 
gorithm for the whole computational process, except for 
the last application of the oracle Uf (compulsory if one 
wishes to distinguish between the two states). It is eas- 
ily seen, because in all the steps of the computation the 
quantum state is an arbitrary superposition of computa- 
tional states of amplitudes +1/ VN and —1/ \/N, without 
any apparent structure, so the probabilities of being in 
|V»o) or \ipo) evolve erratically. The full structure only 
appears when the last oracle Uf is applied, thus provid- 
ing the necessary majorization to be able to distinguish 
the two states with a measurement, but the important 
point is that this is not a majorization arrow, because 
there is no step-by-step majorization. This is stated in 
our second result: 

Result 2: No step-by-step majorization is present 
along the non- efficient parity determination problem. 

Thus, no majorization cycle similar to the one found 
in the preceding section could ever appear. Interestingly 
enough, this is a problem in which quantum computers 
do not provide a better efficiency than classical ones. 



V. ANALYSIS OF THE ADIABATIC 
SEARCHING ALGORITHMS AND OF AN 
ADIABATIC ALGORITHM SOLVING A 2-SAT 
PROBLEM 

We now turn to the analysis of the quantum adiabatic 
searching algorithm, observing the effects of a change of 
path between the initial and the problem hamiltonian 
under the majorization's point of view [l^. We see that 
those paths leading to optimality in the algorithm lead 
as well to step-by-step majorization, while the converse 
is not necessarily true. A different example of adiabatic 
evolution is analyzed in the last point of this section, 
namely, an adiabatic algorithm solving the 2-SAT "ring 
of agrees" problem. 

The adiabatic model of quantum computation was ini- 
tially introduced in Q , and can be briefly summarized as 
follows. We consider a physical system controlled by a a 
time dependent hamiltonian of the form 



H(s(t)) = (1 - s(t))H + s(t)H p , 



(21) 



where H$ and H p are the initial and problem hamilto- 
nian respectively, and s(t) is a time dependent function 
satisfying the boundary conditions s(0) = and s(T) = 1 
for a given T. The desired solution to a given problem 
is encoded in the ground state of H p . The gap between 
the ground and the first excited state of the instanta- 
neous hamiltonian at time t will be called g{t). Let us 
define g m in as the global minimum of g(t) for t in the 



interval [0, T]. If at time T the ground state is given by 
the state \E ;T), the adiabatic theorem states that if we 
prepare the system in its ground state at t = (which is 
assumed to be easy to prepare) and let it evolve under 
this hamiltonian, then 



|<£ ;2W))| 2 > l 



(22) 



provided that 



max 



dH 1: , 
dt 



< e 



(23) 



where -Hi,o is the hamiltonian matrix element between 
the ground and first excited state, e <C 1, and the max- 
imization is taken over the whole time interval [0, T] 
0, |2^. Because the problem hamiltonian encodes the 
solution to the problem in its ground state, we get the 
desired solution with high probability after a time T. 

This quantum adiabatic evolution method has been 
applied to the searching problem ([H El IH 13). Let 
the initial state be \ipo) = ^= Ylx=i I 2 ")' being N the 
number of entries for the searching, and let the initial 
and problem hamiltonian be Hq = I — | i/^o) (V^o | and 
Hp = I — \m)(m\, being |m) the searched state. This 
scheme leads to different results depending on whether 
we apply the adiabatic condition globally (that is, in the 
whole time interval [0,T]) or locally (at each time t). In 
what follows we will consider these two situations without 
entering into precise details of the involved calculations. 
For further information, we refer the reader to j2|| [2tJ 
and references therein. 



A. Analysis of the fastest global adiabatic evolution 

Let us suppose that we demand the adiabatic condi- 
tion given in equation (|23fl to be satisfied globally in the 
whole interval [0, T]. This does not involve any particu- 
lar restriction on the form of s(t), so we can then choose 
s(t) = t/T, leading to a linear evolution of the hamilto- 
nian. Under these circumstances it can be proven |27| 
that the global adiabatic condition is verified provided 
that 



TV 

T> — 

e 



(24) 



Hence, the algorithm needs O(N) time to hit the solu- 
tion with appreciable probability, so the global adiabatic 
searching does not lead to an increasing efficiency with re- 
spect to a classical searching (in contrast with the square 
root speed-up of Grover's quantum algorithm). 

Let us call P+(t) the probability of being at the 
searched state at time t and similarly P-{t) the prob- 
ability of being at any different state from the desired 
one at time t. Note that, because of the symmetry of 
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the problem, P_ (t) will be exactly the same for all those 
quantum states differing from the solution state all along 
the computational process. In order to analyze majoriza- 
tion, we recall the set of inequalities given in equation 
to be satisfied at each majorizing time step. It is easy 
to see that the maximum probability at all times will 
be P+(t), while the rest of the probabilities will remain 
smaller to this quantity all along the computation and 
equal to P-(t). In order to gain simplicity we have an- 
alyzed in detail the behaviour of the two non-trivial cu- 
mulants P + (t) and P + (t) + P_ (t) , as the rest of them will 
not lead to different conclusions from the ones emerging 
from our study. 

We have performed numerical simulations in the fastest 
allowed case (T = ^) and have found the time evolution 
for the two cumulants. The results for e = 0.2 and TV = 
32 are shown in Fig. [5J 
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FIG. 2: Global adiabatic evolution with parameters e = 0.2, 
N = 32 and T = 160. The solid line corresponds to the time 
evolution of P+{t) and the dashed line that of P+(i) + P-(t). 

From our numerical analysis we arrive at our third re- 
sult: 

Result 3: A naive adiabatic quantum searching pro- 
cess does not produce an optimal algorithm neither veri- 
fies step-by-step majorization. 

This property is clearly seen as the two cumulants de- 
crease in time for some time steps, thus not verifying 
majorization. 



B. Analysis of the local adiabatic evolution 

The preceding global adiabatic method can be im- 
proved if we apply the adiabatic condition given in equa- 
tion (|23|1 locally. That is, let us divide the interval [0, T] 
into many subintervals and let us apply ll'.'iH to each of 
the subintervals individually. Taking the limit of the size 



of the subintervals going to zero, we find that the adia- 
batic restriction has to be fulfilled locally at each time 
t: 

W- e v< - (25) 

This is a less demanding condition than l|23|l (if ED is sat- 
isfied, so is 1251 but the inverse is not necessarily true), 
and means that the adiabatic evolution must be satisfied 
at each infinitesimal time interval. It can be shown (see, 
for example, |27j 1 that proceeding in this way the func- 
tion s(t) must have a precise form which is given by the 
relation 

t = ^- = ('arctan(v / A ir - l(2s - 1)) + arctan(v / A f - 1)) 

(26) 

We can observe this dynamic evolution from Fig. |3 in 
the case of e = 0.2 and N = 32. The local adiabatic pro- 
cess implies that the smaller the energy gap between the 
ground and first excited states, the slower the evolution 
of the hamiltonian (and viceversa). 




5 10 15 20 25 30 35 40 
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FIG. 3: Local adiabatic evolution. The driving hamiltonian, 
with e = 0.2 and N = 32. 

With this information it can be proven as well j2]j 
that the evolution time for the algorithm to succeed with 
appreciable probability is, in the limit where JV»1. 

T = ^-VN . (27) 

Hence, in the case of local adiabatic evolution the com- 
putational process takes 0(v~N) time, just as in Grover's 
quantum searching algorithm, obtaining an square root 
speed-up with respect to a classical searching. 

Defining again P + (t) and P_ it) in the same way as in 
Sec. V.A, we can restrict ourselves to the analysis of the 
two non-trivial cumulants P+(t) and P+(t) +P_(t) in or- 
der to observe the evolution of majorization in the quan- 
tum process. We have numerically solved the equations 
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for e = 0.2 and N = 32, and have found the evolution of 
the two quantities, which is given in Fig. 0] 
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FIG. 4: Local adiabatic evolution with parameters e = 0.2, 
N = 32 and T — 44. The solid line corresponds to the time 
evolution of P+(i) and the dashed line that of P+(t) + P-(t). 

From the numerical analysis the following result 
emerges: 

Result 4: A local adiabatic searching algorithm is op- 
timal and verifies step-by-step majorization. 

This result can be straightforwardly verified since the 
set of inequalities of J5J) are satisfied step-by-step, ac- 
cording to the plot in Fig. 0] Due to this behavior, the 
whole computational process might lead to a majoriza- 
tion cycle, such as the one observed in Sec. Ill, as long 
as the previous preparation of the initial quantum state 
of the computation at time t = leads to a stcp-by-step 
reverse majorization. This turns to be always possible, 
for example, by applying a set of Hadamard gates to the 
log AT-qubit quantum state |0 . . . 0) (assume that N is a 
power of 2, for simplicity), which would efficiently pre- 
pare a superposition of all the possible quantum states 
together linked to a reverse majorization arrow. Fur- 
thermore, this quantum process leads to an increasing 
efficiency with respect to a classical searching, exactly in 
the same fashion as Grover's original quantum searching 
algorithm pj. 



C. Analysis of slower global adiabatic evolutions 

Let us now address the situation of global adiabatic 
evolutions which are not necessarily tight in time, that is, 
extremely slow time variations for the hamiltonian, much 
slower than the minimum necessary for the adiabatic the- 
orem to hold. In the case we are dealing with, that im- 
plies the consideration of the case in which T > —, i.e., 
the adiabatic inequality (|24fl is not tight. This case is 



not very relevant from a computational point of view be- 
cause the hitting time is not the minimum possible, but 
we think it is worthwhile to be studied also from the point 
of view of majorization theory in order to have a more 
complete picture of how majorization really works in this 
kind of quantum algorithms. 

We have performed again numerical analysis for the 
time evolution of the two non-trivial cumulants, for e = 
0.2, N = 32, and T = 320, 480 (in both cases bigger than 
f = 160). The results are plotted in Fig. 0and Fig. © 
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FIG. 5: Global adiabatic evolution with e = 0.2, N = 32, and 
T = 320. The solid line corresponds to the time evolution of 
P+(t) and the dashed line that of P+(t) + P-(t). 
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FIG. 6: Global adiabatic evolution with e = 0.2, N = 32, and 
T = 480. The solid line corresponds to the time evolution of 
P+(t) and the dashed line that of P+(t) + P-(t). 

From these two plots, we observe with Fig. [21 that a 
step-by-step majorization appears as long as the evolu- 
tion of the hamiltonian becomes slower and slower. Phys- 
ically, this means that the probability of "jumping" to 
the first excited state diminishes as long as the evolution 
is performed with a very small velocity, thus satisfying 
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better the assumptions of the adiabatic theorem. This 
leads to processes in which there is no quantum speed- 
up but there is indeed a majorization arrow. However, 
these processes are not optimal in time, as we can always 
find faster quantum algorithms for solving the problem. 
Consequently, we arrive at the following conclusion: 

Result 5: Non-optimal quantum algorithms may 
present step-by-step majorization. In particular, step-by- 
step majorization may appear in global adiabatic search- 
ing processes for a slow enough evolution rate. 

It follows that step-by-step majorization cannot be a 
sufficient condition for quantum speed-up. 

We can get some further intuition of the set of results 
presented in this section. Adiabatic quantum searching 
algorithms can be understood (in the limit of large N) 
as a rotation from the initial state to the marked state as 
long as the adiabaticity (either global or local) of the evo- 
lution is conveniently satisfied (see j2?| for details). The 
difference between the global and local conditions turns 
out to be the evolution rate of the rotation angle: local 
adiabatic evolution imposes a rotation at constant rate 
(as in the original Grover's algorithm) whereas global 
adiabatic evolution does not. Because of this rotational 
picture, step-by-step majorization is verified as long as 
the quantum state remains in the instantaneous ground 
state all along the computation. We can now understand 
our results in a finer way. Global adiabatic evolution is 
not a strong enough condition for adiabaticity, thus we 
only see step-by-step majorization when the evolution 
is really slow, in which case the quantum state adiabati- 
cally rotates towards the solution because it remains very 
close to the ground state of the instantaneous hamilto- 
nian. Local adiabatic evolution is a stronger condition 
for adiabaticity the quantum state remains always very 
close to the instantaneous ground state, thus performing 
the rotation towards the solution which in turn involves 
step-by-step majorization. 



D. A further example: a 2-SAT quantum adiabatic 
algorithm solving the "ring of agrees" problem 

Let us consider now a different example of adiabatic 
quantum computation, namely, an adiabatic quantum al- 
gorithm solving the 2-SAT "ring of agrees" problem, as 
stated in 0. As long as 2-SAT can be efficiently solved 
by a classical algorithm in a time 0(poly(n)) (being n 
the number of bits) |30l |. quantum computation can do 
no better than classical computation in this case. The 
problem hamiltonian H p is now a sum of hamiltonians 
involving each of the different clauses of the 2-SAT prob- 
lem, whereas the initial hamiltonian Hq is such that its 
ground state is again an equal superposition of all the 
possible states of the computational basis. The "ring of 
agrees" problem over n bits is defined such that clause j 
acts on bits j and j + 1 where j runs from 1 to n and bit 



n + 1 is identified with bit 1 . Each clause is an "agree" 
clause, which means that 00 and 11 are the satisfying 
assignments. The eigenvectors of the hamiltonian associ- 
ated with clause j are the computational states, in such 
a way that those which "agree" in qubits j and j + 1 
have zero energy (ground states) whereas those which 
"disagree" have energy one. Because the problem hamil- 
tonian is a sum of the n hamiltonians for the n "agree" 
clauses, its ground states are |0)|0) ■ • • |0), |1)|1) • ■ ■ |1), or 
any linear combination of them. 

We have made an analysis of the adiabatic quantum 
algorithm in which the interpolation between the ini- 
tial and problem hamiltonian is linear, s(t) = t/T, for 
the case of 4 qubits (Hibert space of dimension 16), and 
choosing T = 10. The evolution of the 15 cumulants is 
plotted in Fig. □ 




FIG. 7: Evolution of the 15 cumulants in the "ring of agrees" 
for a Hilbert space of dimension 16 



Note the similarity between Fig. Hand Fig. El in Sec. 
V.A, namely, no step-by-step majorization is present in 
the evolution because of the oscilatory behaviour of the 
cumulants. Both plots represent quantum algorithms 
which do not improve classical computation and which 
probabilities share the same behaviour under the point 
of view of majorization. Our observation is then the fol- 
lowing: 



Result 6: A quantum adiabatic algorithm solving the 
2-SAT "ring of agrees" problem does not improve clas- 
sical computation, neither verifies step-by-step majoriza- 
tion. 



This result reinforces the ones already found with re- 
spect to adiabatic scarchig algorithms. 
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VI. ANALYSIS OF A QUANTUM WALK IN 
CONTINUOUS TIME WITH EXPONENTIAL 
ALGORITHMIC SPEED-UP 

The extension of classical random walks to the quan- 
tum world has been widely studied, yielding two dif- 
ferent models of quantum random walks, namely, those 
which operate in discrete time by means of using a "coin 
operator" 0, HI] and those based on a hamilto- 
nian evolution in continuous time 0, 13, HH ■ Regarding 
the discrete time model of quantum random walk two 
interesting algorithmic results have been found so far, 
namely, an exponentially faster hitting time in the hyper- 
cube with respect to the classical random walk and 
a quantum searching algorithm achieving the Grover's 
quadratic speed-up [33 • The first of these examples does 
not provide any algorithmic speed-up, as there exists a 
classical algorithm that solves the hitting problem in the 
hypcrcube exponentially faster than the naive classical 
random walk, that is, in a time 0(poly(logiVY) where 
N is the number of nodes of the graph (see On 
the other hand, the second of these examples shows al- 
gorithmic advantage with respect to any possible clas- 
sical computation. The analysis of the quantum ran- 
dom walk searching algorithm shows that the quantum 
evolution can be understood as an (approximate) rota- 
tion of the quantum state in a two-dimensional Hilbcrt 
space which is exact in the limit of a very large database 
(see [37j for details), resembling the original proposal of 
Grover's searching algorithm which can be decomposed 
exactly in a two-dimensional Hilbert space. This rota- 
tional structure of the evolutin implies step-by-step ma- 
jorization when approaching the marked state, exactly 
in the same fashion than the usual Grover's searching 
algorithm (already analyzed in [l6j'). 

In this section we restrict ourselves to the continu- 
ous time model of quantum walk and analyze a recently 
proposed quantum algorithm based on a quantum walk 
on continuous time solving a classically hard problem 
£| . Here we restrict ourselves to briefly sketch the main 
points and ideas of both the problem setting and its quan- 
tum efficient solution, since the whole development of the 
algorithm is not the purpose of the present paper. We 
address the interested reader to 01. 



A. Setting of the problem 

The problem we wish to solve is defined through a 
graph built in the following way (see 0): suppose we are 
given two balanced binary trees of height n with the 2™ 
leaves of the left tree identified with the 2™ leaves of the 
right tree in a simple way, as shown in Fig. |SJ A way 
of modifying such a graph is by connecting the leaves 
by a random cycle that alternates between the leaves of 
the two trees, instead of identifying them directly. An 
example of such a graph is provided in Fig. 

Suppose that the edges of such a graph are assigned 




OUT 



FIG. 8: A possible graph constructed from two binary trees 
with n = 3. 




OUT 



FIG. 9: An alternative graph constructed from two binary 
trees with n — 3 (connection between the leaves is made 
through a random cycle). 

a consistent coloring (that is, not two edges incident in 
the same vertex have the same color), and that the ver- 
tices are each one given a different name (with a 2n-bit 
string, so there are more possible names than the ones 
assigned). We now define a black box that takes two in- 
puts, a name a given as a 2n-bit string and a color c, 
and acts in the following way: if the input name a corre- 
sponds to a vertex that is incident with an edge of color 
c, then the output corresponds to the name of the vertex 
joined by that edge; if a is not the name of a vertex or 
a is the name of a vertex but there is no incident edge 
of color c, the output is the special 2n-bit string 11 ... 1, 
which is not the name of any vertex. 

Now, the problem we wish to solve reads as follows: 

Given a black box for a graph such as the one previously 
described, and given the name of the IN vertex, find the 
name of the OUT vertex. 
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In |j] it was proven that no classical algorithm can 
transverse a graph such as the one in Fig. [§] in poly- 
nomial time, given such a black box. Furthermore, an 
explicit contraction of a quantum algorithm based on a 
continuous time quantum walk on the graph that suc- 
ceeds in finding the solution for this oracular problem in 
polynomial time was given. 

B. Quantum algorithm 

The quantum algorithm of can be briefly summa- 
rized as follows: consider the (2n + 2)-dimensional sub- 
space spanned by the states 



achieved with a very high probability. In [4| the authors 
prove that the succeeding time is polynomial in n. 

C. Analysis of the quantum algorithm 

We have numerically simulated this quantum walk for 
the particular case of n = 4, and have plotted the time 
evolution of the probability of success in Fig. El We 
observe that the numerical result fits with the prediction 
that the time the algorithm takes in achieving the OUT 
node is polynomial. 

1 i 1 1 1 1 1 1 



I col j) 



E 



(28) 



3 aGcolumn j 



where Nj = 2 j if < j < n and Nj = 2 2n+1 ~i if 
n + 1 < j < 2n + I. We will call this subspace the "col- 
umn subspace" , and each state of the basis is an equally 
weighted sum of the states corresponding to the vertices 
lying on each column of the graph. We now define a 
hamiltonian acting on this subspace in the following way: 



0.6 



0.4 



0.2 



/ , -i tti w ,m f 1 < j < n — 1 , n + 1 < j < 2n 
(colj|//|col(j + l)) = { ^ - J n ~ ~ J - 

(29) 

with hermiticity of H giving the other nonzero matrix 
elements. The action of this hamiltonian in the graph is 
nothing but promoting transitions between adjoint ver- 
tices, so a quantum walk on the graph (on the whole 
Hilbert space) generated by this hamiltonian is equiva- 
lent to a quantum walk on the line (on the column sub- 
space). Consequently, from now on we only focus our 
attention in the quantum walk in the line generated by 
the hamiltonian of l|29|) . Moreover, it can be proven that 
given the structure of the graph in the form of a black 
box such as the one already described, our hamiltonian 
can be efficiently simulated 

The quantum walk works as follows: at first the "wave 
packet" will be precisely localized at the IN vertex (the 
initial state will be |col 0)). Due to unitary time evo- 
lution, it will initially spread out through the different 
vertices at the left hand side of the graph (those be- 
longing to the left binary tree), but after a short time 
(once half the graph has been transversed) it will begin 
to spread through the vertices on the right hand side, 
interfering constructively in the OUT vertex as the time 
goes on. Physically, this is nothing but a wave propaga- 
tion. Should we wait more time, the wave packet would 
come back to the entrance, and the process would be sim- 
ilarly repeated again. Actually, due to the "defect" of the 
hamiltonian in the middle vertices, it can be shown that 
the transmission through the central columns is not of 
100 percent (thus providing interferences in long enough 
time scales), but high enough for the OUT node to be 
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FIG. 10: Quantum walk algorithm. Probability of finding the 
exit, for n = 4. 

It is easy to observe that, in order to analyze majoriza- 
tion, for the case of n = 4 there are only 10 non-trivial 
probabilities to be studied. This is so due to the fact that 
all the states of the whole Hilbert space belonging to the 
same column always share the same probability ampli- 
tude. The relevant quantities to be studied are then the 
probabilities of being at each column state normalized 
to the number of nodes belonging to that column, that 
is, the probability of being in one node of each column. 
There are then 2n + 2 different probabilities to be con- 
sidered at each time step. Given only this 10 quantities, 
we are able to calculate the whole set of 62 cumulants 
corresponding to all the sums of sorted probabilities, ac- 
cording with equation J5J). In order to make the figures 
as clear as possible we only plot 10 of these quantities 
in Fig. 1111 corresponding to the cumulants arising from 
the sorted probabilities when only one node per column 
is considered. The rest of the cumulants can be shown 
to have a similar behavior to that of the ones appearing 
in Fig. El 

We have also numerically simulated the algorithm in 
the case of a bigger graph, namely, in the case n = 10. 
In this case there are 2rt + 2 = 22 different probabilities 
to be considered at each time step. Proceeding in the 
same way than in the case n = 4 (that is, not plotting 
all the cumulants, but the only the sorted sum of these 
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FIG. 11: Quantum walk algorithm. Time evolution of the 
ten cumulants when one node per column is considered, for 
n = 4. 



22 probabilities), we obtain a similar behaviour as in the 
case for n = 4, as is shown in Fig. IT21 



E 




FIG. 12: Quantum walk algorithm. Time evolution of the 
22 cumulants when one node per column is considered, for 
n = 10. Note that most of the quantities apparently collapse 
when making the plot because of the small difference in the 
probabilities given the big size of the graph. 

Therefore, we arrive at the following conclusion: 

Result 7: The continuous time quantum walk 
transversing a classically hard graph follows a step-by- 
step majorization cycle all along the computation until it 
reaches the OUT node. 

It is worth remarking as well that the time the algo- 
rithm spends reversely majorizing the probability distri- 
bution is about half of the time of the whole computation. 
The physical reason for this behavior is clear, as this is 
the time the "wave packet" spends spreading over the 



binary tree on the left hand side, thus leading to a de- 
structive interference part. It is note worthy that such a 
destructive interference indeed strictly follows a step-by- 
step reverse majorization of probabilities (satisfying the 
inequalities given in equation © for the case of reverse 
majorization). Furthermore, we see combining Fig. I1UI 
and Fig. ^Jthat the growing of the probability of success 
is linked to a step-by-step majorization. Physically, this 
is the part in which the algorithm constructively inter- 
feres into the OUT node once the wave packet is approx- 
imately in the right hand side binary tree. We see that 
this constructive interference follows a majorization ar- 
row, thus verifying step-by-step the inequalities given in 
equation J2J. Actually, the observed majorization cycle 
reminds us the one already found in the quantum algo- 
rithm of Sec. Ill, but in this case we have numerically 
checked that the present cycle does not seem to follow the 
rules of natural majorization. Complcmcntarily, we have 
also observed that the probability amplitudes follow the 
interesting rule that those belonging to even columns are 
real and those belonging to odd columns arc imaginary. 

The deterministic search by quantum random walk 
heavily exploits the column structure of the problem. 
The register works on a superposition of columns, that is 
of states belonging to the same column with equal weight. 
It is then natural to ask whether a step-by-step majoriza- 
tion cycle operates also at the level of columns. The idea 
behind this analysis corresponds to accept that the final 
measurement will filter each columns as a whole. The 
result of the measurement would correspond to deter- 
mining a particular column. The subtle point here is to 
find to what extend the success of finding the OUT state 
is related to the column structure of the algorithm. We 
have numerically considered the column amplitudes for 
n = 4 and n = 10 with a total of 9 and 21 cumulants to 
be calculated respectively from the sorted probabilities 
at each time step of being at each column of the graph. 
In Fig. 1131 and Fig. ^]we plot the result, which shows 
that there does not exist a majorization cycle when the 
final measurement is made on columns. 

The conclusion is that deterministic quantum walks 
cleverly exploit the column subspace structure of the 
problem to achieve step-by-step majorization on the in- 
dividual states. 



VII. CONCLUSION: A MAJORIZATION 
PRINCIPLE 

We can now collect all the results found in the anal- 
ysis of majorization in the quantum algorithms we have 
studied so far and synthesize an emerging principle un- 
derlying all of them. There are a total of nine empirical 
observations about step-by-step majorization: 

• Presence of step-by-step majorization in Grover's 
quantum searching algorithm |l6[ . 
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FIG. 13: Quantum walk algorithm. Time evolution of the 
nine cumulants when the column measurement is considered, 
for n = 4. 
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FIG. 14: Quantum walk algorithm. Time evolution of the 21 
cumulants when the column measurement is considered, for 
n = 10. 



• Presence of a natural step-by-step majorization cy- 
cle in quantum phase estimation algorithms [l7|. 

• Presence of natural step-by-step majorization cycle 
in the quantum algorithm for finding hidden affinc 
functions (Sec. III.C). 

• Absence of step-by-step majorization in an optimal, 
yet non-efficient, quantum algorithm for solving the 
parity problem (Sec. IV. C). 

• Absence of stcp-by-step majorization in naive 
global quantum adiabatic searching algorithms 
(Sec. V.A). 

• Presence of step-by-step majorization in appropri- 
ate local quantum adiabatic searching algorithms 
(Sec. V.B). 



• Emergence of step-by-step majorization for a slow 
enough evolution rate in quantum adiabatic search- 
ing algorithms (Sec. V.C). 

• Absence os step-by-step majorization in an adia- 
batic quantum algorithm solving the 2-SAT "ring 
of agrees" problem (Sec. V.D). 

• Presence of step-by-step majorization cycle in a de- 
terministic quantum walk on a graph solving a clas- 
sically hard problem (Sec. VI. C). 

Note that our results concerning the analysis of adia- 
batic quantum algorithms can have an alternative valid 
interpretation: according to Fig. [21 and Fig. [7| we see 
that the part of the processes which docs not obey step- 
by-step majorization only occurs when the probability of 
succes has almost achieved its highest value. We are then 
led to the consideration that absence of majorization only 
appears once the algorithms have already constructed the 
right solution, having already done their job. A redefi- 
nition of the algorithm by stopping the process once the 
probability of the winner is maximum would lead us to 
affirm that step-by-stcp majorization is naturally present 
in quantum algorithms by adiabatic evolution according 
to the evidence presented here. This new interpretation 
does not alter our final result. 

Adiabatic algorithmic processes do lead as well to a 
reverse majorization of the probability distribution in or- 
der to efficiently prepare the initial quantum state of the 
computation. This can be efficiently performed by a set 
of Hadamard gates (which produce step-by-step natu- 
ral reverse majorization). Nevertheless, this remark does 
not only hold for the adiabatic paradigm. The usual 
formulation of Grover's algorithm in terms of quantum 
gates needs as well of a preparation of the initial quan- 
tum state which can be carried out exactly in the same 
way. Similarly, all algorithms accomodate to a reverse 
majorization-majorization cycle. We could argue that 
the initial stcp-by-step reverse majorization procedure 
at the beginning of the quantum algorithms is some- 
how trivial, as it only involves the application of (for 
instance) a set of hadamard gates (with the exception of 
the quantum walk algorithm, in which the initial step- 
by-step reverse majorization is by no means trivial in the 
sense we state here as it is carried by the structure of the 
graph). The fact that the quantum evolution accomo- 
dates to a step-by-step reverse majorization-majorization 
cycle is reminiscent of the reversibility of these quantum 
algorithms. 

All the results found so far suggest that a stcp-by- 
step reverse majorization-majorization cycle seems to be 
a necessary condition for efficiency in quantum compu- 
tational processes, although not sufficient. This can be 
promoted to a principle: 

Majorization Principle: 

Optimal quantum algorithms must follow a majoriza- 
tion cycle. 
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This principle fits well with all the observed results 
given so far. Note that those processes which are not 
optimal do not necessarily follow the majorization cycle 
pattern: the case of the optimal algorithm solving the 
parity problem does not, while the extremely slow and in- 
efficient but majorizing adiabatic processes do. Step-by- 
step majorization may be viewed as a strong irreversibil- 
ity condition for success probability necessary for optimal 
quantum algorithms. 

All our results are also consistent with a stronger 
statement, namely that both step-by-step majorization 
and large entanglement complement each other and are 



needed for exponential speed-up. Entanglement brings 
the genuine quantum mechanical tool which has to be 
used in an optimal way, that is verifying step-by-step 
majorization. 
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